arXiv:1502.04808vl [math.CA] 17 Feb 2015 


New patterns of travelling waves in 
the generalized Fisher-Kolmogorov equation 


Pavel Drabek 1 

Department of Mathematics and 
N.T.I.S. (Center of New Technologies for Information Society) 
University of West Bohemia 
P.O. Box 314 

CZ-306 14 Plzen, Czech Republic 
e-mail: pdrabekOkraa. zcu. cz 


and 

Peter Takac 2 

Institut fur Mathematik 
Universitat Rostock 
UlmenstraBe 69, Haus 3 
D-18055 Rostock, Germany 

e-mail: peter. takac@uni-rostock. de 


February 18, 2015 


Abstract. We prove the existence and uniqueness of a family of travelling waves in a degenerate 
(or singular) quasilinear parabolic problem that may be regarded as a generalization of the semilinear 
Fisher-Kolmogorov-Petrovski-Piscounov equation for the advance of advantageous genes in biol¬ 
ogy. Depending on the relation between the nonlinear diffusion and the nonsmooth reaction function, 
which we quantify precisely, we investigate the shape and asymptotic properties of travelling waves. Our 
method is based on comparison results for semilinear ODEs. 

1 The work of Pavel Drabek was supported in part by the Grant Agency of the Czech Republic 
(GACR) under Grant #13 — 00863S. 

2 The work of Peter Takac was supported in part by Deutsche Forschungsgemeinschaft (DFG, Ger¬ 
many) under Grant # TA 213/16-1. 


1 



Running head: Travelling waves in the FKPP equation 

Keywords: Fisher-Kolmogorov equation, travelling waves, 

nonlinear diffusion, nonsmooth reaction function, 
comparison principle 

2010 Mathematics Subject Classification: Primary 35Q92, 35K92; 

Secondary 35K55, 35K65 


2 


1 Introduction 


The purpose of this article is to investigate a very basic pattern formation in a reaction- 
diffusion model, namely, travelling waves. The model is the favorite Fisher-KPP equation 
(or Fisher-Kolmogorov equation) derived by R. A. Fisher [6] in 1937 and first mathe¬ 
matically analyzed by A. Kolmogorov, I. Petrovski, and N. Piscounov [11] in the 
same year. However, these original works ([6, 11]) consider only linear diffusion and (suffi¬ 
ciently) smooth (nonlinear) reaction. In our present work, we allow for both, a nonlinear 
diffusion operator (with a (p — l)-homogeneous quasilinear elliptic part, 1 < p < oo) and 
a nonsmooth reaction function of Holder class C ,0, “ _1 (R) with 1 < a < 2, 


( 1 . 1 ) 


du 

~dt 


d_ 

dx 



du 

dx 



f(u ), (i,i)elx R + , 


u(x, t) = U(x — ct ) 


for some speed c G R . 


Here, R + = [0,oo), 1 < p < oo, d : R — I K is a positive continuous function, and 
/ : R —y R is a continuous, but not necessarily smooth function of a “generalized” Fisher- 
KPP type (specified below). For p = 2, this equation describes a deterministic version of 
a stochastic model for the spatial spread of a favored gene in a population, suggested by 
R. A. Fisher [6]. 


Closely related situations and solution ideas have already been explored, e.g., in P. 
Drabek, R. F. Manasevich, and P. Takac [2] and in Y. Sh. Il’yasov and P. 
Takac [10], where existence, uniqueness, and stability of phase transition solutions in 
a Cahn-Hilliard-type model are investigated. More precisely, we study the interaction 
between the (nonlinear) diffusion and the (nonsmooth) reaction; in paticular, we explore 
their influence on the formation and the shape of a travelling wave connecting two stable 
(spatially constant) steady states. Our main result, Theorem 3.1, contains the existence 
and uniqueness of such a family of travelling waves (parametrized by a spatial shift). 
This result is crucial for establishing long-time front propagation (convergence) in non¬ 
linear parabolic equations towards a travelling wave; see, e.g., D. G. Aronson and 
H. F. Weinberger [1], P. C. Fife and J. B. McLeod [5], and F. Hamel and N. 
NADIRASHVILI [7] for p = 2 (the original semilinear Fisher-KPP equation), E. FEIREISL 
et AL. [4] for any p G (1, oo), and Q. Yl and J.-N. ZHAO [16] for p G (2, oo) only. 

In order to prove Theorem 3.1 we use a phase plane transformation (cf. J. D. Mur¬ 
ray [13], §13.2, pp. 440-441) to investigate a nonlinear, first order ordinary differential 
equation with an unknown parameter c G R, see (3.6). Since this differential equation 
is supplemented by homogeneous Dirichlet boundary conditions at both end-points, this 
boundary value problem is over determined. We find a unique value of c for which this 
problem has a unique positive solution. The nonlinearity in this differential equation does 
not satisfy a local Lipschitz condition, so, due to the lack of uniqueness of a solution, 
the classical shooting method cannot be applied directly. The novelty of our approach is 
to overcome this difficulty: We take advantage of monotonicity properties of / and treat 
the differential equation in problem (3.6) as an initial (at —1) or terminal (at +1) value 
problem. We thus derive various comparison principles that compensate for the lack of 
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uniqueness in the shooting method. Finally, we obtain new shapes for travelling waves 
by asymptotic analysis near the end-points =pl. 

Since the seminal paper by A. M. Turing [14] has appeared in 1952, the (linear) 
diffusion has been known to have a destabilizing effect on stable steady states in a reaction- 
-diffusion equation with a smooth reaction function. This effect leads to the formation 
of new, more complicated patterns, for instance, in morphology used in Turing’s work 
[14]. Among other things, we will demonstrate that the analogues of Turing’s findings 
do not apply universally to our nonlinear problem setting. (Our situation is similar to, 
but not identical with that treated in [14].) We will determine a simple, exact relation 
between the constants p and a (1 < p < oo and 1 < a < 2), when a destabilizing effect 
occurs (1 < a < p < 2) and when it does not occur (1 < p < a < oo). Loosely speaking, 
this effect depends on the product 

(“diffusion”) x (“smoothness”) = - • a < (>) 1. 

p 

Although our results bear resemblance to Turing’s observations [14] made in a branching 
(bifurcation) setting, his destabilizing effect occurs precisely for p = a = 2. In this case, 
we do not have any branching phenomenon in our model in the sense that the values of 
our travelling wave stay in the open interval (—1,1). In contrast, we speak of a branching 
phenomenon exactly when the travelling wave attains one of the values =pl at a (finite) 
spatial point x T i G R. 

In his monograph [12], §14.9, pp. 424-430, J. D. Murray studies “nonexistence 
of stable spatial patterns for scalar equations in one dimension with zero flux bound¬ 
ary conditions” . In particular, the following interesting conclusion is derived on p. 426: 
Large diffusion prevents spatial patterning in reaction diffusion mechanisms with zero flux 
boundary conditions. Drawing an analogue of this conclusion to our setting, we will show 
(Theorem 3.1) that this is the case if p < a. To be more precise, we investigate mono¬ 
tone travelling waves in the degenerate (or singular) second-order parabolic problem of 
a “generalized” Fisher-KPP type (1.1). For p = 2 and / smooth, extensive studies of 
travelling waves can be found in P. C. Fife and J. B. McLeod [5] and J. D. Murray 
[ 13], §13.2, pp. 439-444. 

Recall that d : R —> R is continuous and positive. The function / : R —y R is 
assumed to be continuous, such that /(±1) = /(s 0 ) = 0 for some —1 < s 0 < 1, together 
with f(s) > 0 for every s G (—1, So), /(s) < 0 for every s G (so, 1), and 

/ r 

d(s) 1/(p_1) f(s) ds > 0 whenever — 1 < r < 1. 

Typical examples for the diffusion coefficient d(s) > 0 are 

(dl) d(s) = 1 for all s G R, which yields the p-Laplacian in eq. (1.1), 1 < p < oo. 

(d2) p — 2 and d(s) = (p'(s) > 0 for all s G R, which yields the porous medium operator 
u i —y -7^z<p(u) on the right-hand side of eq. (1.1), where the function p : R —> R is 
assumed to be monotone increasing and continuous. 
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An important special case of the reaction function / is /(s) = F'(s) where F : R —» R 
is the “generalized” double-well potential F(s) = F a (s) = ^|s 2 — 1|" for s G R, where 
a G (1, oo); hence, so = 0 and 

(1.3) f(s) = f a (s) = |s 2 — l|“ _2 (s 2 — l)s for s G R. 

Notice that, if 1 < a < 2 then the function /(s) is only (a — l)-Holder continuous at 
the points s — ±1, but certainly not Lipschitz continuous. If the diffusion coefficient 
d : R — y R (already assumed to be a positive continuous function) is also even about 
zero, that is, d(—s) = d(s) for all set, then (1.3) forces 

(1.4) G(l) = J d(s) 1/(p-1) f(s) ds = 0 


in formula (1.2). 

This article is organized as follows. Our main results are collected in Section 3. In the 
next section (Section 2), the travelling wave problem for the Fisher-KPP equation (1.1), 
the quasilinear ODE (2.1), is transformed into an equivalent problem for an expression 
y : (—1,1) —y R, the semilinear ODE (2.5) with the boundary conditions (2.7), cf. prob¬ 
lem (3.6). The substituted unknown expression y is a simple function of the travelling 
wave U(x) and its derivative U'(x), thus yielding a simple differential equation for the 
travelling wave U, on one hand. On the other hand, in eq. (2.5), the unknown function 
y depends solely on U G (—1,1) as an independent variable, i.e., y = y(U). This prob¬ 
lem, with a monotone, non-Lipschitzian nonlinearity, is solved gradually in Sections 4 
(existence) and 5 (uniqueness), respectively. Finally, an important special case (with a 
Fisher-KPP-type reaction function), which is a slight generalization of (1.3), is treated in 
Section 6. 


2 Preliminaries 


Assuming that the travelling wave takes the form u(x,t) = U(x — ct ), (x,t) G R x R + , 
with U : R —> R being strictly monotone decreasing and continuously differentiable with 
U' < 0 on R, below, we are able to find a first integral for the second-order equation for 
U: 


( 2 . 1 ) 



d{U) 


dU 


dx 



+ C ^ _/ ( f/ ) = °’ Xe 


Following the standard idea of phase plane transformation (U, V) for the p-Laplacian (cf. 
[3, Sect. 1]), we make the substitution 


V d = - d(U) 


dU 


dx 


V -2 


dU 

dx 


> 0 , 
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whence 

( 2 . 2 ) 


d U _ ( V \ 1/(p ~ 1} 

dx ~~ \d(U)) 


and consequently look for V = V(U) as a function of U G (—1,1) that satisfies the 
following differential equation obtained from eq. ( 2 . 1 ): 


dV d U , dU n 

n 77 'T^ + c n- f(U) = 0, 16 

d U dx ax 


that is, 


(2.3) 


dV ( V V /(p - 1} ( V V /(p - 1} 

dc/'teJ ~ c \m> 


f(U) = o, U g (-1,l). 


Finally, we multiply the last equation by d(U) 1 ^ p make the substitution y = f V p ' > 0, 
where p 1 = p/(p — 1) G (1, oo), and write r in place of U, thus arriving at 


V ? - c y 1/P - d (r ) 1/(p ~ 1] Hr) = 0 , r G (-1,1). 

p' dr 


This means that the unknown function y : (—1,1) —> (0, oo) of r, 


(2.4) 


y = V p ' = d(UY 


dU 

dx 


p 

> 0 , 


must satisfy the following differential equation: 

(2.5) ^ =P {c(y + ) 1/p + g(r)) , r G (-1,1), 

where y + = ma x{y, 0 } and g(r) = f d(r) 1 ^ p ~ 1 ' 1 f(r) satisfies the same hypotheses as /: 


g : M —* R is a continuous, but not necessarily smooth function, such that g(±l) = 
g(so) = 0 for some —1 < s 0 < 1, together with g(s) > 0 for every s G (—l,s 0 ), 
g(s) < 0 for every s G (so, 1), and (1.2), i.e., 

/ r 

g(s) ds > 0 whenever — 1 < r < 1 . 


Since we require that U = U (x) be sufficiently smooth, at least continuously differentiable, 
with U'(x) —* 0 as x —> ±oo, the function y = y{r) must satisfy the boundary conditions 

(2.7) 2 /( — 1) = 2/(1) = 0 . 


The following remark on the value of (7(1) ( > 0) is in order. 
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Remark 2.1 Since the integrand s t-G g(s) '= d(s) 1 ^ p ~ 1 ' > f(s) : (—1,1) —> R in the 
function G(r), dehned in (1.2) for r E (—1,1), is continuous and absolutely integrable 
over (—1,1), we conclude that G : [—1,1] —» R is absolutely continuous. In particular, 
ineq. (2.6) forces G(l) > 0. We will see later that the case G(l) = 0 guarantees the 
existence of a stationary solution to problem (1.1), i.e., c = 0, whereas the case G(l) > 0 
renders a travelling wave , i.e., c ^ 0; more precisely, with c < 0, cf. Theorem 3.1 below. 
Indeed, both, the stationary solution (for c = 0) and the travelling wave (for c < 0) will 
be obtained from eq. (2.5) by means of the transformation defined by eqs. (2.2) and (2.4). 


3 Main Results 


We assume that d : R —> R is a positive continuous function. In fact, we need only 
d : [—1,1] —> R to be continuous and positive. Let us recall that the function g(r) = 
d(r) 1 /( p_1 ) f(r) of r E [—1,1] is assumed to satisfy the hypotheses formulated in the 
previous section (Section 2) before Remark 2.1, in particular, ineq. (2.6). 

Our main result concerning the travelling waves u(x,t) = U(x — ct), (x,t) G R x R + , 
c G R, for problem (1.1) is as follows. 


Theorem 3.1 Let G(l) > 0. Then there exists a unique number c* G R such that 
problem (1.1) with c = c* possesses a travelling wave solution u(x,t) = U(x — c*t), 
(x,t) G RxR + , where U = U{f) is a monotone decreasing and continuously differentiable 
function on R taking values in [—1,1]. Furthermore, we have c* < 0 and the set {(gR: 
U'(f) < 0} is a nonempty open interval (xi,x_i) C R with 

(3.1) lim [/(£) = 1 and lim !/(£) = — 1. 

£-Kaa)+ £-Kz-i)- 


Proof. The existence and uniqueness of c* G R follow from the transformation of 
eq. (2.1) into eq. (2.5) with the boundary conditions (2.7) and subsequent application of 
Theorem 3.4 below to the boundary value problem (2.5), (2.7). We get also c* < 0. 

Inserting y = y c * : (—1,1) —> R into eq. (2.4) we obtain 

J V(r) — 2/ c *(r) 1 ^ / = y c *('r’Y P ~ 1 ^ p >0 for r G (—1,1); 

\v(-l) = V(l) = 0. 


We insert V = V{U) as a function of U G (—1,1) into eq. (2.2), thus arriving at 


dC7 
dx 


fV(U )\ 1/(p “ 1} 

{W)) 


< 0 


for U G (-1,1). 
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Separation of variables above yields 


dx= ~ {v$) ' {r ' iU for u e (-1. i) 


Finally, we integrate the last equation to arrive at 




(3.2) 


x{U) = x(0) — 


cl(r) 

W: 


v(p-i) 


dr for U G (—1,1), 


where x(0) = Xq G R is an arbitrary constant. We remark that V (r) —y 0 as r —y ±1. 
Consequently, both monotone limits below exist, 


(3.3) x_i = lim x(U) and aq = lim x(U) , 

{/->•(-1)+ c/-K+i)- 

and satisfy —oo < aq < x 0 < x_i < +oo. Thus, the function x : (—1,1) —?• (aq,a;_i) C R 
is a diffeomorphism of the open interval (—1,1) onto (aq, x_i) satisfying ^ < 0 in (—1,1) 
together with (3.3). This implies the remaining part of our theorem, especially the limits 
in (3.1). ■ 


ft depends on the asymptotic behavior of the function g = g(r), r G (—1,1), near the 
points ±1 whether U' < 0 holds on the entire real line R = (— oo, +oo) or else U' < 0 
on a nonempty open interval (aq,a;_i) C R and U' = 0 on its nonempty complement 
R\(aq, a;_i) with aq > —oo and/or x_i < +oo. More precisely, we assume that g : R —* R 
has the following asymptotic behavior near ±1: 

There are constants G (0, oo) such that 


(3.4) 


sM 


lim , 

r—>1— (1 — r)' 


= -7o 


and 


g(r) 


lim 

r — 1+ (1 + r) 


— 7o 


For instance, if g(s) = d(s ) 1 ^ p f(s) for s G (—1,1) is as in Remark 2.1 and f(s) = 

f a (s ) is dehned by (1.3), then we have a = 1 + n / ± . 

We have the following conclusions for the limits x^i = lim t /_ l .( =Fl ) ± x(U), where x = 
x{U) is a solution of 


(3.5) 


dx 

dU 


d{U) \ 1/(p_1) 
V(U)J 


< 0 


for —1<[/<1, 


cf. eq. (2.2), given explicitely by formula (3.3). Notice that, in this notation, x : (—1,1) —> 
(aq,x_i) is the inverse function of the restriction of U to the interval (aq,x_i) in which 
U' < 0, i.e., x{r) = U~ l {r ) for —1 < r < 1. 


Theorem 3.2 Assume that the limits in eq. (3.4) hold. Then, 







(i) in case 1 < p < 2 we have x\ = —oo (x_i = +oo) if and only if 7 + > p — 1 
( 7 ~ > p — 1 , respectively); whereas 

(ii) in case 2 < p < oo we have x\ = — oo (x_i = +oo) f/ 7 + > p — 1 ( 7 “ > p — 1 , 
respectively). 


This theorem is an easy combination of Theorem 3.1 above with Corollary 6.3 and 
Remark 6.4 in Section 6 . 

The following corollary to Theorem 3.2 for the linear diffusion case (p = 2) is obvious: 


Corollary 3.3 Let p = 2. Then x\ (x_i) is finite if 7 + < 1 (7 < 1, respectively) and 
x\ = —00 (x_i = + 00 ) if 7 + > 1 ( 7 “ > 1, respectively). 


The results stated above hinge upon the existence and uniqueness results for eq. (2.5) 
with the Dirichlet boundary conditions (2.7), i.e., for the following boundary value prob¬ 
lem: 


^7 = p'{c{y + ) 1/p + g(r)) , re (- 1 , 1 ); 

?/(—!) — 2 /( 1 ) — 0 , 


with the parameter c 6 R to be determined. We recall that eq. (2.5) has been obtained 
from eq. (2.2) by means of the substitution in (2.4). 

The following result for problem (3.6) is of independent interest (cf. [4, 16]). 


Theorem 3.4 Let G(l) > 0. Then there exists a unique number c* G M such that 
problem (3.6) with c = c* has a unique solution y = y c * : (—1,1) —$■ R. Furthermore, we 
have c* < 0 and y c * > 0 in (— 1 , 1 ). 


Proof. The existence and uniqueness of c* G R follow from Corollary 5.5 (with the 
existence established before in Corollary 4.7). We get also c* < 0. ■ 


4 Existence Result for c* 


In order to verify that only the case c < 0 can yield a travelling wave for problem (1.1), 
we prove the following simple lemma for the boundary value problem (3.6). 
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Lemma 4.1 Let c G R and yo(r) = p' G(r) for —1 < r < 1. 

(i) Let c = 0. Then the function y 0 (r ) = p' G{r) of r G [—1,1] is a solution to 
problem (3.6) if and only if G( 1) = 0 . 

(ii) U c > 0 (c < 0, respectively ) then every solution y : [—1, 1] —> R to the initial 
value problem for eq. (2.5) with the initial condition y(— 1) = 0 satisfies y > yo (y < yo) 
throughout (— 1 , 1 ]. 


Proof. Part (i) is a trivial consequence of (2.6) combined with the properties of g. 

Part (ii): First, let c > 0. Then, clearly, y' > y' 0 = p' g throughout (—1,1), together 
with y' ^ y' 0 in (—l,r) for every r G (—1,1), which forces y > yo throughout (—1,1], by 
a simple integration of eq. (2.5) over the interval [—l,r]. 

The case c < 0 is analogous, by reversing the (nonstrict) inequalities. ■ 


For c > 0, Lemma 4.1 has the following obvious consequence. 


Corollary 4.2 If either (i) c = 0 and G(l) > 0 or else (ii) c > 0 and G(l) > 0, 
then eq. (2.5) possesses no solution y : [—1,1] —* R satisfying the boundary conditions 
y(—1 ) = 0 and y( 1) < 0 . 


We conclude that it suffices to investigate the case c < 0 for finding a travelling wave 
to problem ( 1 . 1 ). 


Remark 4.3 Given any c < 0, the nonlinearity y H > c{y + ) l ^ p : R — y R is nonincreasing 
and, thus, satisfies a one-sided Lipschitz condition, which guarantees that eq. (2.5) with 
the initial condition y(— 1 ) = 0 possesses a unique solution y : [— 1 , 1 ] —> R, by two well- 
known existence and uniqueness theorems from Ph. Hartman’s monograph [8], Theorem 
2.1, p. 10 (existence due to Peano), and Theorem 6.2, p. 34 (one-sided uniqueness), 
respectively. This solution depends continuously on c < 0, by [ 8 , Theorem 2.1, p. 94], 
y = y C: and it may change sign. Furthermore, for every fixed r G (—1,1], the function 
c i —y ?/ c (r) : M_ —y R is monotone increasing, by [8, Corollary 4.2, p. 27]. Here, we have 

denoted R_ (—oo,0]. We investigate this dependence after this remark. 

Even if we do not need to exploit the case c > 0 any more, we still need to consider 
the terminal value problem for eq. (2.5) with the terminal condition y( 1) = 0, i.e., the 
“backward” initial value problem for eq. (2.5) with y( 1) = 0. Unfortunately, if c < 0, 
we do not have a uniqueness result of the kind described above for this problem. To 
overcome this difficulty, in the next section (Section 5) we establish a comparison result 
in Proposition 5.4(b). 
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More detailed results concerning existence, uniqueness, monotone dependence on cer¬ 
tain parameters, and other qualitative properties of solutions to the initial value problem 
for eq. (2.5) with the initial condition y(— 1) = 0 will be established in Corollary 5.3. 

In view of Remark 2.1 and Corollary 4.2 above, from now on we assume 

(4.1) G(l) = fd(s)^ns)as>0 

in formula (1.2). Equivalently, we have i/o(l) = p' G{ 1) > 0 
a unique constant c < 0 such that the (unique) solution y c 
value problem for eq. (2.5) with the initial condition y c {— 1) - 
condition y c iX) = 0. To find this constant, let us define 

(4.2) c* d = inf (cei: y c ( 1) > 0} . 

As expected, we will show that precisely c* is the desired value of the constant c; see 
Corollaries 4.7 and 5.5. 

We need the following two technical lemmas. 


We will show that there is 
[— 1 , 1 ] —» R to the initial 
0 satisfies also the terminal 


Lemma 4.4 Let c G M be arbitrary. 

(i) There is some 5 = S(c) G (0, 2) such that y c (r ) > 0 holds for all r G (—1, —1 + 5). 

(ii) If z/ c (l) > 0 then we have also y c {r) > 0 for all r G (—1,1). 


Proof. Part (i): On the contrary, if no such number 6 G (0,2) exists, then there is a 
sequence {r n }™ =1 C (—1,1) such that r„ \ -1 as n / oo, and y c (r n ) < 0 for all n G N. 
Given any fixed n G N, we cannot have y c (s) < 0 for all s G (—l,r n ] since, otherwise, by 
eq. (2.5), the function y c would satisfy 

(4.3) ^=p'g(r), r G (—1, r n ), 

dr 

that is, y c = V o > 0 throughout [—1, r n ], by a simple integration of eq. (4.3) over the inter¬ 
val [—1,7’]. This conclusion would then contradict our choice y c (r n ) < 0. Consequently, 
there is a point r' n G (—l,r„) such that ?/ c (r(J > 0. ft follows that there is another point 
r" G (r' n , r n \ such that y c (r") < 0 and y' c (r ") < 0. Since r n \ -1 as n / oo, we have also 
r" —* — 1 as n —> oo. But now, substituting r = r" in eq. (2.5), we arrive at g[r") < 0 
for every n G N. A contradiction with our hypothesis on g , namely, g(s) > 0 for every 
s G (—1, s 0 ), is obtained for every n G N sufficiently large, such that r" G (—1, s 0 ). 

We have proved the existence of the number 6. 

Part (ii): On the contrary, let us assume that the set {r G (—1,1) : ydj") < 0} is 
nonempty. We denote by r G (—1,1) the smallest number r G (—1,1) such that y c (r) = 0. 
Hence, — 1 < — l + 5<r<l and y c (r) > 0 holds for every r G (—1, r). We have y c (r) = 0 
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and, therefore, p'g(r) = y'Jfr) < 0, by eq. (2.5). Our hypothesis on the nodal point s 0 
of the function g thus forces r > so- Hence, we must have y' c (r) = p'g(r) < 0 for all 
r G (r, 1) C (so, 1). We conclude that y c (r) < 0 holds for all r G (r, 1], which contradicts 
our hypothesis y c ( 1) > 0. 

We have proved that y c {r) > 0 holds for all r G (—1,1) as desired. ■ 


Lemma 4.5 Assume that c < 0 and G(l) > 0. Let y : [—1,1] —y R be a solution 
to eq. (2.5) satisfying y(— 1) = y(r 0 ) = 0 for some r 0 G (—1,1) and y(r ) > 0 for all 
r G (—l,r 0 ). Then we have ro > So and y(r) < 0 holds for all r G (r 0 ,1]. 


Proof. On the contrary, suppose that —1 < r 0 < s 0 . Then y'(r 0 ) = p' g(r 0 ) > 0 
holds by our hypotheses on g. However, this is contradicts our hypothesis on y requiring 
y{r) > 0 = y(r 0 ) for all r G (— l,r 0 ) and, thus, forcing y'(r 0 ) < 0. This proves r 0 > s 0 . 

Consequently, we have g(r) < 0 for all r G (ro, 1) which entails y'(r ) = p' g(r) < 0 for 
all r G (r 0 , 1). Integration over the interval [r 0 ,r] yields the desired result, i.e., y(r) < 0 
holds for all r G (r 0 , 1], ■ 

Proposition 4.6 We have —oo < c* < 0 and y c *iX) = 0. 


Proof. First, let us recall that the terminal value y c ( 1) G M is a monotone increasing, 
continuous function c i—z/ c (l) : —y M of the parameter c G M_. These properties yield 

immediately — oo < c* < 0 and, if c* > — oo then also y c *{X) = 0. 

Thus, it remains to prove c* > — oo. By contradiction, let us suppose that there is 
a sequence of numbers c n G (—oo, 0), such that y Cn ( 1) > 0 for every n — 1, 2, 3, ... and 
c n \ —oo as n oo. Then ycnir) > o must hold for all r G (—1,1], by Lemma 4.4, 
Part (ii). 

Next, let z n : [so, 1] —> M be the (unique) solution of the initial value problem 
a 7 

(4.4) = p'c n {zf) 1/p , r G (s 0 , 1); z n (s 0 ) = y Cn {s 0 ) ■ 

The monotone (increasing) dependence of the solution y c : [—1,1] —y M on the right-hand 
side of eq. (2.5) (essentially due to E. Kamke; see, e.g., M. W. Hirsch [9, p. 425] or W. 
Walter [15, Chapt. Ill, §10], Comparison Theorem on p. 112) guarantees the following 
comparison result: 0 < 7/ Cri (r) < z n (r) for all r G [so, 1]. However, separating variables in 
eq. (4.4) and integrating, we obtain 

(■ z n (r )) 1/p ' = (^(s 0 )) 1/p ' + c n (r - s 0 ) for all r G [s 0 , 1] . 

Consequently, by the comparison result, we have also 

0 < (yc n (r)) 1/p ' < (; y Cn (s 0 )) 1/p ' + c n (r - s 0 ) for all r G [s 0 ,1] . 
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Setting r = 1 and recalling y Cri ( 1) > 0, we observe that {y Cn {sfi)) l ^ p ' > —c n ( 1 — s 0 ) holds 
for every n = 1,2,3,.... Since also y Cn < y 0 holds throughout the interval [—1,1], by 
Lemma 4.1, Part (ii), we conclude that — c n (l — So) < (2/o(so)) 1,/p , a contradiction to 
— c n /* oo as n / A oo. We have proved c* > —oo as desired. ■ 

Lemma 4.5 and Proposition 4.6 entail the following obvious corollary. 

Corollary 4.7 There is some number c 6 R, say, c = c* (< 0), such that problem (3.6) 
possesses a solution y : [—1,1] —>■ M. This solution is unique, given by y = y c *, and it 
satisfies y(r) > 0 for all r e (—1,1). 


For this particular speed c = c*, the uniqueness of y c * has been discussed in Re¬ 
mark 4.3. The sole existence of y c * can be found also in R. Engukja, A. Gavioli, and 
L. Sanchez [3, Theorem 4.2, p. 182] with a sketched proof. 


5 Uniqueness Result for c* 


For the sake of reader’s convenience, in this section we establish a few comparison results 
(weak and strong) which provide important technical tools. These results are essentially 
due to E. Kamke, as we have already mentioned in the proof of Proposition 4.6 above. 
They will enable us to establish the uniqueness of c* specified in Corollary 4.7 (see Corol¬ 
lary 5.5 below). 


Definition Let us consider the differential equation (2.5) on some nonempty open in¬ 
terval (a,b) C (—1,1) C M, regardless an initial condition at a or a terminal condition 
at b. 

We say that a continuously differentiable function y : [a, b] —> M is a subsolution to 
equation (2.5) on (a, b ) if and only if 

(5.1) < p' {c(y + ) 1/p + g(r)) , re(a,b). 

Similarly, we say that a continuously differentiable function y : [a, b] —* M is a supersolu¬ 
tion to equation (2.5) on (a, b ) if and only if 

(5.2) ^ >p'{c(y + ) 1/p + g(r)) , r e (a,b). 


Remark 5.1 We remark that the terminology for sub- and supersolutions introduced 
above fits the case when eq. (2.5) is supplemented by an initial condition for y(a) at the 
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initial point r = a. However, if eq. (2.5) is supplemented by a terminal condition for 
y(b) at the terminal point r = b, it would be logically more appropriate to switch the 
inequalities (5.1) and (5.2) in the definitions of sub- and supersolutions; compare Parts 
(a) and (b) of Proposition 5.2 below. Namely, a natural way to view a terminal value 
problem is to revert the (“time”) variable r i —y —r which entails also the reversion of 
inequalities (5.1) and (5.2). However, we will not perform this reversion to avoid possible 
confusion. 


We begin with the following comparison result. 


Proposition 5.2 Let c6i, — 1 < a < b < 1, and assume that y,y : [a, b] —y R is a pair 
of continuously differentiable functions, such that y (y) is a subsolution ( supersolution, 
respectively ) to equation (2.5) on ( a,b ). 

(a) If c < 0 and y(a) < y(a) then we have also y(r) < y(r) for all r £ [a, b\. 

(b) If c > 0 and y(b ) > y(b) then we have also y{r ) > y(r) for all r £ [a, b\. 


Recall our abbreviations M + 


= [0, oo) and 


def 


(-oo,0]. 


Proof. We prove only Part (a); Part (b) is proved analogously. 

Hence, assume that c < 0. We subtract ineq. (5.2) from (5.1), thus arriving at 

(v - V) < p'c (( y + ) 1/p - (; y + ) 1/p ) , r e (a, b). 

Now let us multiply this inequality by (y — y) + , thus obtaining 


(5.3) (y-y) + -^;(y-y)<p'c({y + ) 1/p -{y + ) 1/p ){y-y) + < 0, re(a,b), 

thanks to c < 0 combined with the montonicity of the functions £ H> and £ i —> (£ + ) 1//p 

from M to R + . The last inequality entails 


1 A 

2 dr 


{(y-y) + ) 2 < 0 


for almost every r £ (a, b ), 


which shows that w = ((y — y) + Y is a monotone nonincreasing, nonnegative function on 
[a, b\. Since w(a) =0 by our hypothesis, it follows that w(r) — 0 holds for every r £ [a, b\. 
This completes our proof of Part (a). ■ 


The following corollary on the uniqueness and the monotone dependence on the pa¬ 
rameter c £ M_ of the solution y = y c : [—1,1] —y M to the initial value problem for 
eq. (2.5) with the initial condition y(— 1) = 0 is an easy direct consequence of Proposi¬ 
tion 5.2. 
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Corollary 5.3 (i) Given any c G M_, the initial value problem for eq. (2.5) with the 
initial condition y(— 1) = 0 possesses a unique (continuously differentiable ) solution y = 
Vc ■ [ — 1 ; 1 ] 

(ii) V — oo < ci < C2 < 0 then y Cl < y C2 holds pointwise throughout [—1,1], In 
particular, the function c K > y c { 1) : R_ —> K is monotone increasing. 


Proof. Part (i) (the uniqueness claim) follows trivially from Proposition 5.2, Part (a). 

Part (ii): Let — oo < Ci < C2 < 0 and choose any constant c G [ci, C2]; hence, 
c < 0. Then y Cl ,y C2 '■ [ — 1,1] —> M is a pair of sub- and supersolutions to equation (2.5) 
on (—1,1), respectively, with c chosen above (ci < c < C2) and the initial condition 
y Cl (— 1) = Uc 2 {~ 1 ) = 0. The pointwise ordering y Cl < y C2 throughout [—1,1] now follows 
from Proposition 5.2, Part (a), again. ■ 


We complement Proposition 5.2 with another comparison result. 

Proposition 5.4 Let c G M, a G [so, 1), and assume that y,y : [a, 1] —> M is a pair 
of continuously differentiable functions, such that y (y) is a subsolution ( supersolution, 
respectively) to equation (2.5) on (a, 1). Moreover, assume that 

y(r) > 0 and y(r) > 0 for all r G [a, 1), together with y( 1) = 0 < y( 1). 

Then we have also 

(a) y{r) < y(r) for all r G [a, 1]. 

(b) More precisely, if y( 1) = 0 < y( 1) then y(r) < y{r) holds for all r G [a, 1]. If 

y(l) = y( 1) = 0 then there is some number a' G [a, 1] such that y(r) < y(r ) for all 

r G [a, a') and y(r) = y(r) for all r G [a', 1]. 

The comparison result in Part (b) takes care of the lack of uniqueness in the terminal 
value problem for eq. (2.5) with the terminal condition y( 1) = 0 and c < 0, as announced 
in Remark 4.3. 

Proof of Proposition 5.4. Part (a): We begin by substituting the variable t = 1 — r 
with 0 <f<l — a (< 1 — So) and the functions 

z{t) = ( 2/(1 - t)) 1/p ' and z{t) = (y(l - t)) 1/p ' 
with the derivatives (recall that ^ ^ = 1 ) 

^ z(t) = “ ^ (vi 1 ~ t ))~ 1/p ■ ^y{l ~ t) and 

^7 z{t) = - — (y(l - t))~ 1/p ■ y{l - t) for 0 < t < 1 - a . 

dt p at 
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Hence, inequalities (5.1) and (5.2), respectively, are equivalent with 


(5.4) 

dz 

— > 
d t 


t G (0,1 — a 

and 




(5.5) 

dz 

— < 

dt 


f G (0,1 — a 


We subtract ineq. (5.4) from (5.5), thus arriving at 
(5,6) 4 (- 2 _ 2 ) < - g(l -1) [(z(t))-‘«D] , t £ (0, 1 - o]. 

Now let us multiply ineq. (5.6) by (z — z) + , thus obtaining 


(5,7) (7 - z) + 4 (z - z) < - 9 (1 - f) [(7(i))- 1/(p - 1) - (sM)- 1 ^ 1 ’] (z - z)+ < 0 


for all i e (0,1 
entails 


a], thanks to g(l — t ) < 0 for 0 < t < 1 — a < 1 — Sq. The last inequality 


1 d 

2 di 


(( 7-‘) + ) 2 <0 


for almost every r 6 (0,1 — a), 


which shows that w = f ((z — z) + ) 2 is a monotone nonincreasing, nonnegative function on 
[0,1 — a}. Since w(0) = 0 by our hypothesis, it follows that w(t) = 0 holds for every 
t £ [0,1 — a]. Equivalently, y{r) < y{r ) holds for all r 6 [a, 1], which proves Part (a). 


Part (b): We return to ineq. (5.6). The difference on the right-hand side takes the 
form 


0 < (z(t)) 


rt+W-Vip- 1 ) 


m) 


- 1 / 0 - 1 ) _ 


(5.8) 


p — 1 


[(1 - 6)z(t) + 9z{t)\ p d 9 (. z(t) - z(t )) 


< L ( z(t ) — z(t)) for all t E [1 — b, 1 — a ], 


where b G (a, 1) is arbitrary and the Lipschitz constant L = L(b) G (0, oo) depends on b. 
We may take 

1 / \ ~ p ' 

L = -- min z(t) G (0, oo), 

p — 1 ytG[l—6,1—a] J 

by the inequality 


/ [{1 - 9)z(t) + 9z(t)}~ p ' dd< (z{t))~ p \ f G (0,1 — a]. 

Jo 

Consequently, we can estimate the right-hand side of ineq. (5.6) as follows: 


(5.9) 


-(z-z)<L (. z{t) - z(t )) , t G [1 — b, 1 — a ], 
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where 


L = L(b) = f L ■ sup \g\ , L G (0, oo) 

[- 1 , 1 ] 


e Lt ( z(t ) — z(£)) < 0 for every t G [1 — b, 1 — a] . 


Ineq. (5.9) is equivalent with 

< 5 - 10 > s 

By integration over any compact interval [ti, t 2 ] C [1 — b, 1 — a], we get 
(5.11) «(t 2 ) - z(t 2 ) < e _i(t2_tl) (z(ti) - z(ti)) < 0 . 


First, assume that ]/(l) = 0 < y( 1), that is, z(0) = 0 < z(0). Then, by continuity, also 
z(ti) < z(t\) holds provided t\ G (0, 1 — a) is chosen to be small enough. From ineq. (5.11) 
we deduce that also z(t) < z(t) for all t G [ti, 1 — a]. We have proved z(t) < z(t) for every 
t G [0,1 — a]. The desired inequality y{r) < y{r) for every r G [a, 1] follows immediately. 

Now assume y( 1) = y( 1) = 0, that is, z(0) = z(0) = 0. Let t' G [0,1—a] be the greatest 
number t G [0,1 — a], such that z{t) — z(t) = 0 holds for every t G [0,£']. Equivalently, 
y{r ) — y(r) = 0 holds for every r G [a 7 ,1], where a' = 1 — t' G [a, 1]. If t' — 1 — a then we 
are done. Thus, assume 0 < t' < 1 — a. Then there is a sequence {r n C (t\ 1 — a), 
such that r n \ t' as n /*■ oo and z(r n ) — z(r n ) < 0 for every n — 1,2,3,.... Again, 
from ineq. (5.11) we deduce that also z(t) < z(t ) for all t G [r n , 1 — a]. We have proved 
z(t) < z(t) for every t G (t f , 1 — a]. The desired inequality y(r) < y{r ) for every r G [a, a') 
follows immediately. 

The proof is complete. ■ 


Corollary 4.7, Corollary 5.3, and Proposition 5.4 entail the following important result. 

Corollary 5.5 There is a unique number c G M such that problem (3.6) possesses a 
solution y : [—1,1] —$■ M, namely, c = c* (< 0) defined in eq. (4.2). This solution is given 
by y = y c * and satisfies y(r) > 0 for all r G (—1,1). 


Proof. The existence of c* (< 0) has been established in Corollary 4.7. On the 
contrary to uniqueness, suppose that there is another number c, say, c = c' G M, c' c *, 
such that y = y c > satisfies the boundary conditions y(— 1) = y( 1) = 0. We have y c '{r) > 0 
and y c *(r) > 0 for all r G (—1,1), by Lemma 4.4, Part (ii). By our definition of c* in 
eq. (4.2), we must have c' < c* (< 0). Fixing any constant c G [c', c*], we observe that y c > 
(; y c *, respectively) is a subsolution (supersolution) to equation (2.5) on (—1,1) (with this 
c). 


Corollary 5.3, Part (ii), implies y c ' < y c * throughout [—1,1]. In contrast, Propo¬ 
sition 5.4, Part (b), forces y c * < y c ' on [s 0 , 1 ]. We conclude that y c i = y c * on [s 0 , 1 ]. 
Consequently, both, y c i and y c *, verify eq. (2.5) on the interval (so, 1 ) for both values of c, 
c = d and c = c*. Since also y c >(r) = y c *{r) > 0 for every r G [s 0 ,1), eq. (2.5) with c — d 
and c = c* forces c' = c*. We have proved also the uniqueness. ■ 
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6 Asymptotic Behavior 


In this section we investigate the asymptotic behavior of the (unique) solution y = y c * : 
[—1,1] —> M to problem (3.6) with c = c*. Of course, we assume G(l) > 0 throughout the 
present section. We recall that the existence and uniqueness of both, c* (< 0) and y c *, 
have been obtained in Corollaries 4.7 and 5.5. 

In fact, the asymptotic estimates for r G (—1,1) near the endpoints ±1 obtained in 
this section remain valid in the following more general situation with an arbitary value of 
the parameter c < 0 near the left endpoint —1: We assume that y = y c : [—1,1] —> R is 
the (unique) solution to eq. (2.5) satisfying the initial condition y(— 1) = 0. 

We assume that g{r ) = d(r) 1 ^ p ~ 1 ' 1 f(r ) satisfies also the following hypothesis: 


There are constants 7, y 0 G (0, 00) such that 


( 6 . 1 ) 


lim 

r—>—1 


g(r) 

(1 + r) 7 


= 7o • 


Equivalently, we have 

(6.1') g(r) = (70 + + r ) 7 for r G [-1,1] , 

where r) : [— 1 , 1 ] —)■ M is a continuous function with 7 ( 0 ) = 0 . 

We abbreviate the differential operator 

( 6 . 2 ) (Ay)(r)^‘^--p l ( C (y*ylr + g(r)) , re (-1,1), 

for any function y G C 1 ([—1,1]). We are interested in constructing sub- and supersolutions 
to equation (2.5) as defined in Section 5. Let k G (0, oo) be a constant (to be determined 
later) and let us consider the (nonnegative) function 

(6.3) w K (r) = f n (1 + r ) 1+7 for r G [—1,1] . 

Then we have 

(Aw K )(r) = - p' ( c ( w K ) 1/p + g(r )) 

(6-4) = + y)(l + r ) 7 — p' (cK 1,/p (l + r)^ 1+7 ' )//p + (y 0 + p(r))(l + r) 7 ) 

= [fi)(l + 7 ) ~ P ( 7 o + v( r ))] (1 + r ) 7 — p'cK 1/,p (l + r) <yl+1 ^ p 

for r G (—1,1). Letting r G -1 we observe that the first expression, with the power 
( 1 +r) 7 , dominates the second one, with the power (l+r) (7+7 - )/,p , if anc [ only if 7 < (l+ 7 )/p 
holds, that is, 7 < p'/p = l/(p — 1). Consequently, if this is the case, then, for every 
r G (—1,1) sufficiently close to —1 we have: 
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(Aw K )(r) < 0 provided n > 0 is small enough, and 
(Aw K )(r) > 0 provided n > 0 is large enough. 

We can state these simple facts about w K being a sub- or supersolution to problem (2.5) 
on some interval (— 1 , — 1 + g) as follows (0 < g < 2 ): 

Lemma 6.1 Let c E R and assume that g : [—1,1] —» R satisfies (6.1) with 0 < 7 < 
l/(jp — 1). Then there exist numbers g E (0, 2) and 0 < k < k < 00 , such that 

(i) if 0 < k < k then ( Aw K )(r ) < 0 holds for all r E (—1, —1 + 0 ); and 

(ii) if k < k < 00 then ( Aw K )(r ) > 0 holds for all r E (—1, — 1 + g). 


We remark that all three numbers g, k, and k may depend on c, in general. 

Now, let us fix c = c* given by eq. (4.2) and recall Corollary 5.5. We combine 
Proposition 5.2, Part (a), with Lemma 6.1 to conclude that 

(6.5) w 5 (r) = k (1 + r ) 1+7 < y c *(r) < w K (r) — k(1 + r ) 1+7 

for all r E (-1,-1 + g). We recall that V(U) = [y c *{T)} l ^ p with U — r, by eq. (2.4). 
Consequently, inequalities (6.5) above read 

( 6 . 6 ) (k) 1/p ' (1 + U) {1+l)/p ' < V{U) < (1 + f/) (1+7)/p ' 


for all U E (-1, 
(3.5) 


— 1 + g). I 11 eq. (2.2), that is 
dx _ fd(l 7 )\ 1/(p - 1} 


< 0 , 




for the unknown function x = x(U), we need the following equivalent form of ( 6 . 6 ), 

(6.6') ( k) 1/p (1 + [/)( 1+7 )/p < I/(f/) 1/(p - 1) < ( k) 1/p (1 + U)^ ),p 


for all U E (—1, —1 + g). 

It remains to investigate the case 7 > (1 + 7 )/p, that is, 7 > p'/p = 1 /(p — 1), in 
which the second expression in eq. (6.4), with the power (1 + r) ( ' 1+7 ' )//p , dominates the first 
one, with the power (1 + r) 7 . Here, we assume c < 0. Then, given any k E (0,oo), we 
have (thanks to c < 0 ) 


( Aw K )(r) > 0 for every r E (—1,1) sufficiently close to —1. 


We can state this simple fact about w K being a supersolution to problem (2.5) on some 
interval (— 1 , — 1 + g) as follows (0 < g < 2 ): 
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Lemma 6.2 Let c < 0 and assume that g : [—1,1] —y R satisfies (6.1) with 7 > l/(p — 1). 
Then, given any k e (0, 00 ), there exists a number g e (0, 2), such that 

(. Aw K )(r ) > 0 holds for all r 6 (— 1 , — 1 + g). 


We remark that, again, the number g may depend on k and c, in general. 


Now, let us fix c = c* given by eq. (4.2) and recall Corollary 5.5. We combine 
Proposition 5.2, Part (a), with Lemma 6.1 to conclude that 

(6.7) Vc*(r) < w K (r) = k (1 + r ) 1+7 

for all r 6 (- 1,-1 + g). We recall that V(U) = [y c *(r)] 1/,p with U — r, by eq. (2.4). 
Consequently, ineq. (6.7) above reads 

(6.8) V(U) < k 1/p ' (1 + U) (1+7)/p ' 


for all U e (—1, —1 + g). In eq. (3.5), that is 

dx _ fd(U) \ 1/(p_1) 

w ~ \nuj) 


< 0 , 


-1 < U < 1, 


for the unknown function x = x(U), we need the following equivalent form of ( 6 . 8 ), 
( 6 . 8 ') V{U) 1/[ - p - l) < k 1/p (1 + U) ( 1 + T/p 


for all U G (—1, —1 + g). 


Corollary 6.3 Let c = c* (< 0) and assume that g : [—1,1] —y R satisfies (6.1). Then 

we have the following conclusions for the limit x_i l = lim^.i x(U) ( < + 00 ) of the 
solution x = x{U) to the differential equation (3.5): 

(i) In case 1 < p < 2, the limit x_i is finite if and only if 0 < 7 < p—1. ( Equivalently, 
X-i = +00 if and only if 7 > p — 1 .) 

(ii) In case 2 < p < 00 , the limit x_i = +00 if 7 > p — 1. 

Consequently, we have x_i = +00 whenever 0 < p — 1<7<oo. 


The proof of Part (i) (Part (ii), respectively) follows directly from eq. (3.5) combined 
with ineq. ( 6 . 6 )' (ineq. ( 6 . 8 )'). 

In the linear diffusion case, i.e., for p = 2, x_i is finite if 7 < 1, whereas x_i = +00 if 
7 > 1 . 

Remark 6.4 Under a condition similar to (6.1), but imposed at the right endpoint as 

def 

r —y +1, we obtain conclusions analogous to those in Corollary 6.3 for the limit x\ = 
limjy-n-i x(U) ( > - 00 ). 
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